Scalar Quarkonium Masses and Mixing with the Lightest Scalar 

Glueball 



W. LeeE and D. Weingarten 
IBM Research, P.O. Box 218, Yorktown Heights, NY 10598 

Abstract 

00 
On 

£>Y We evaluate the continuum limit of the valence (quenched) approximation 

to the mass of the lightest scalar quarkonium state, for a range of different 
quark masses, and to the mixing energy between these states and the lightest 
scalar glueball. Our results support the interpretation of /o(1710) as com- 
posed mainly of the lightest scalar glueball. 
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Evidence that /o(1710) is composed mainly of the lightest scalar glueball is now given 
by two different sets of numerical determinations of QCD predictions using the theory's 
lattice formulation in the valence (quenched) approximation. A calculation on GF11 [[Ij] of 
the width for the lightest scalar glueball to decay to all possible pseudoscalar pairs, on a 
16 3 x 24 lattice with j3 of 5.7, corresponding to a lattice spacing a of 0.140(4) fm, gives 
108(29) MeV. This number combined with any reasonable guesses for the effect of finite 
lattice spacing, finite lattice volume, and the remaining width to multibody states yields a 
total width small enough for the lightest scalar glueball to be seen easily in experiment. For 
the infinite volume continuum limit of the lightest scalar glueball mass, a reanalysis |2j of 
a calculation on GF11 ||, using 25000 to 30000 gauge configurations, gives 1648(58) MeV. 
An independent calculation by the UKQCD-Wuppertal [|J collaboration, using 1000 to 3000 
gauge configurations, when extrapolated to the continuum limit according to Ref. |2],|5| yields 
1568(89) MeV. A more recent calculation using an improved action || gives 1630(100) MeV. 
The three results combined become 1632(49) MeV. A phenomenological model of the glueball 
spectrum which supports this prediction is discussed in Ref. 0. 

Among established resonances with the quantum numbers to be a scalar glueball, all are 
clearly inconsistent with the mass calculations except /o(1710) and /o(1500). Between these 
two, / (1710) is favored by the mass result with largest statistics, by the combined result, 
and by the expectation || that the valence approximation will lead to an underestimate of 
the scalar glueball's mass. Refs. |l|J^| interpret / (1500) as dominantly composed of strange- 
antistrange, ss, scalar quarkonium. A possible objection to this interpretation, however, is 
that / (1500) apparently does not decay mainly to states containing an s and an s quark ||. 
In part for this reason, Ref. [|Hj] interprets /o(1500) as composed mainly of the lightest scalar 
glueball and /o(1710) as largely ss scalar quarkonium. A second objection is that while the 
Hamiltonian of full QCD couples quarkonium and glueballs, so that physical states should be 
linear combinations of both, mixing is not treated quantitatively in Ref. |l| . In the extreme, 
mixing could lead to /o(1710) and /o(1500) each half glueball and half quarkonium. 

Using the valence approximation for a fixed lattice period L of about 1.6 fm and a range 
of different values of quark mass, we have now calculated the continuum limit of the mass 
of the lightest scalar qq states and the continuum limit of the mixing energy between these 
states and the lightest scalar glueball. The continuum values are obtained by extrapolation 
of results obtained from four different values of lattice spacing. For the two largest lattice 
spacings we have also done calculations on lattices with L of about 2.3 fm. Preliminary 
versions of this work are reported in Refs. P, |TT| . |T2"| . 

Our results provide answers to the objections to the interpretation of / (1710) as largely 
the lightest scalar glueball. For the valence approximation to the infinite volume continuum 
limit of the ss scalar mass we find a value significantly below the valence approximation scalar 
glueball mass. This prediction rules out, in our opinion, the possibility of identifying [l(| 
/o(1500) as primarily a glueball and / (1710) as primarily ss quarkonium. Our calculation of 
glueball-quarkonium mixing energy, combined with the simplification of considering mixing 
only among the lightest discrete isosinglet scalar states, then yields a mixed /o(1710) which is 
73.8(9.5)% glueball and a mixed / (1500) which is 98.4(1.4)% quarkonium, mainly ss. The 
glueball amplitude which leaks from / (1710) goes almost entirely to the state / (1390), 
which remains mainly {uu + dd) For (mm + dd) /\/2 in the rest of this article we will use 
the abbreviation nn, normal- ant inormal. We find, in addition, that /o(1500) acquires an nn 
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amplitude with sign opposite to its ss component suppressing, by interferences, the state's 
decay to KK final states. Assuming SU(3) flavor symmetry before mixing for the decay 
couplings of scalar quarkonium to pairs of pseudoscalars, the KK decay rate of / (1500) is 
suppressed by a factor of 0.39(16) in comparison to the rate of an unmixed ss scalar. This 
suppression is consistent, within uncertainties, with the experimentally observed suppression. 

Our calculations, using Wilson fermions and the plaquette action, were done with en- 
sembles of 2749 configurations on a lattice 12 2 x 10 x 24 with (3 of 5.70, 1972 configurations 
on 16 3 x 24 with (3 of 5.70, 2328 configurations on 16 2 x 14 x 20 with (3 of 5.93, 1733 config- 
urations on 24 4 at (3 of 5.93, 1000 configurations on 24 2 x 20 x 32 with (3 of 6.17, and 1003 
configurations on 32 2 x 28 x 40 with (3 of 6.40. For (3 of 5.70, 5.93, 6.17 and 6.40 the corre- 
sponding values of lattice spacing a are, respectively, 0.140(4) fm, 0.0961(25) fm, 0.0694(18) 
fm, and 0.0519(14) fm. The smaller lattices with (3 of 5.70 and 5.93, and the lattices with 
/3 of 6.17 and 6.40 have periods in the two (or three) equal space directions of 1.68(5) fm, 
1.54(4) fm, 1.74(5) fm, 1.66(5) fm, respectively, and thereby permit extrapolations to zero 
lattice spacing with nearly constant physical volume. 

These values of lattice spacing, and conversions from lattice to physical units in the 
remainder of this article, are determined |||!J from the exact solution to the two-loop zero- 
flavor Callan-Symanzik equation for A-^a with of 234.9(6.2) MeV determined from 



MS MS 

. / ( m n 

MS 



the continuum limit of (A^-=a)/(m p a) in Ref. |13 |. For (3 from 5.70 to 6.17, the ratio 



(Aj^a)/(m p a) in Ref. [H| was found to be constant within statistical errors, thus our results 
are, within errors, almost certainly the same as those we would have obtained by converting 
to physical units using values of m p a. We chose to convert using A^-a, however, since 



Ref. jnj did not find m p a at (3 of 6.40, which would be needed for our present calculations. 
For each ensemble of gauge fields, we evaluated correlation functions using random 



sources built from quark and antiquark fields following Ref |jT2|l . Averaged over random 
sources, these correlation functions become 

<?//(*) = £< f(x,t)f -(0,0) >, 

X 

Cgsit) =]T< g(x,t)s(0,0) > . (1) 

x 

Here / is either p, s or g. The quantities p(x,t) and s(x,t) are, respectively, the smeared 
pseudoscalar and scalar operators of Ref. fl2| built from a quark and an antiquark field and 
g(x,t) is the smeared scalar glueball operator of Ref. ||. 

Fitting the the large-t behavior of the diagonal correlators to the asymptotic form 

c ff(t) Z f exp{-m f at), (2) 

we obtained the masses, in lattice units, m p a, m s a, and m g a and field strength renormaliza- 
tion constants Z p , Z s and Z g . From the large-t behavior of the off-diagonal correlator, for 
m s close to m g , 

C gs (t) — ► J~Zg~Z~sEa^2 ex P(~ m g a \t — t'\ — m s a \t'\) (3) 
v 

we then found the glueball-quarkonium mixing energy, in lattice units, Ea. Eqs. (|2|) and 
(0) have been simplified by omitting terms arising from propagation around the lattice's 
periodic time boundary. 
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For the two lattice with (3 of 5.93, Figure [3] shows the scalar quarkonium mass as a 
function of quark mass m q a, defined to be (2k) _1 — (2k c ) -1 . Here k is the hopping constant 
and k c is the critical hopping constant at which the pseudoscalar mass m p goes to zero. We 
determined k c from a fit of m 2 to a quadratic function of 1/n. The solid lines in Figure [l| are 
quadratic fits to the scalar mass as a function of quark mass which we used to interpolate 
to the strange quark mass. The strange quark mass we chose to be the value which yields, 
in physical units, a pseudoscalar mass squared of 2m 2 K — m%, where tuk and m n are the 
observed neutral kaon and pion masses, respectively. As shown by the figure, for the lattice 
16 2 x 14 x 20 with L of 1.54(4) fm the scalar mass as a function of quark mass flattens out 
as quark mass is lowered toward the strange quark mass and then appears to begin to rise 
as the quark mass is decreased still further. This feature is absent from the data at f3 of 5.93 
for the lattice 24 4 with L of 2.31(6) fm and is thus a finite- volume artifact. It is present in 
the data at (3 of 5.70 with L of 1.68(5) fm, at f3 of 6.17 with L of 1.74(5) fm, and at (3 of 
6.40 with L of 1.66(5) fm, but absent in the data at (3 of 5.70 with L of 2.24(7) fm. 

The pseudoscalar mass squared m 2 , for all values of lattice spacing, we found to be 
nearly a linear function of 1/k and nearly independent of lattice period. The difference in 
m p a between the two lattice at (3 of 5.70 and between the two lattice at (3 of 5.93 was in all 
cases less than 0.5%. 

For L near 1.6 fm, Figure |2] shows the ss scalar mass in units of A-^ as a function of 

lattice spacing in units of A linear extrapolation of the mass to zero lattice spacing 

gives 1322(42) MeV, far below our valence approximation infinite volume continuum glueball 
mass of 1648(58) MeV. For the ratio of the ss mass to the infinite volume continuum limit 
of the scalar glueball mass we obtain 0.802(24). Figure |2| shows also values of the ss scalar 
mass at (3 of 5.70 and 5.93 with L of 2.24(7) and 2.31(6) fm, respectively. The ss mass 
with L near 2.3 fm lies below the 1.6 fm result for both values of lattice spacing. Thus 
the infinite volume continuum ss mass should lie below 1322(42) MeV. We believe our data 
rule out the interpretation of /o(1500) as mainly composed of the lightest scalar glueball 
with / (1710) consisting mainly of ss scalar quarkonium. For comparison with our data, 
Figure |2| shows the valence approximation value for the infinite volume continuum limit of 
the scalar glueball mass and the observed value of the mass of /o(1500) and of the mass of 
/o(1710) The uncertainties shown in the observed masses in units of arise mainly from 

the uncertainty in A^-. 

Figure [3] shows the quarkonium-glueball mixing energy as a function of quark mass for 
the two different lattices with (3 of 5.93. For neither lattice does there appear to be any 
sign of the anomalous quark mass dependence found in Figure [1]. The mixing energies at 
different quark masses turn out to be highly correlated and depend quite linearly on quark 
mass. For (3 of 5.70, 6.17 and 6.40 the mixing energy behaves similarly and, in particular, 
also depends quite linearly on quark mass. Thus it appears that the mixing energy can be 
extrapolated reliably down to the normal quark mass m n , defined to be the quark mass at 
which m p becomes m^. At (3 of 5.70, the mixing energy ratio E(m n )/ E(m s ) is 1.222(34) for 
L of 1.68(5) fm and 1.194(45) for L of 2.24(7) fm. For the data at (3 of 5.93, this ratio is 
1.183(32) for L of 1.54(4) fm and 1.153(56) for L of 2.31(6) fm. Thus the ratio has at most 
rather small volume dependence and seems already to be near its infinite volume limit with 
L around 1.6 fm. 

Figure |] shows a linear extrapolation to zero lattice spacing of quarkonium-glueball 
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mixing energy at the strange quark mass E(m s ) and of the ratio E(m n )/E(m s ). The limiting 
value of E(m s ) is 43(31) MeV and of E{m n )/E{m s ) is 1.198(72). 

We now combine our infinite volume continuum value for E(m n )/E(m s ) with a simpli- 
fied treatment of the mixing among valence approximation glueball and quarkonium states 
which arises in full QCD from quark-antiquark annihilation. The simplification we introduce 
is to permit mixing only between the lightest scalar glueball and the lowest lying discrete 
quarkonium states. We ignore mixing between the lightest glueball and excited quarkonium 
states or multiquark continuum states, and we ignore mixing between the lightest quarko- 
nium states and excited glueball states or continuum states containing both quarks and 
glueballs. 

Excited quarkonium and glueball states and states containing both quarks and glueballs 
are expected to be high enough in mass that their effect on the lowest lying states will be 
much smaller than the effect of mixing of the lowest lying states with each other. On the other 
hand, the shift in the glueball mass arising from mixing with discrete states is, according 
to the systematic version of the valence approximation described in Ref. mainly a 



one-quark-loop correction to the valence approximation while mixing with continuum states 
is mainly a multiquark loop correction. For low energy QCD properties it is expected 
multiquark loop corrections will be significantly smaller than single loop corrections. 

The structure of the Hamiltonian coupling together the scalar glueball, the scalar ss and 
the scalar nn isosinglet becomes 

E(m s ) V2rE( 
E{m„) m S o 



V2rE(m s ) rn 



nn ■ 



Here r is the ratio E{m n ) / E{m s ) which we found to be 1.198(72), and m g , m S s and m n n 
are, respectively, the glueball mass, the ss quarkonium mass and the nn quarkonium mass 
before mixing. 

The three unmixed mass parameters we will take as unknowns. We will also take E{m s ) 
as an unknown since the fractional error bar on our measured value is large. These four 
unknowns can now be determined from four observed masses. To leading order in the 
valence approximation, with valence quark-antiquark annihilation turned off, corresponding 
isotriplet and isosinglet states composed of u and d quarks will be degenerate. For the 
scalar meson multiplet, the isotriplet (uu — dd) / \/2 state has a mass reported by the Crystal 
Barrel collaboration to be 1470(25) MeV f|. Thus we take m nTi to be 1470(25) MeV. In 
addition, the Crystal Barrel collaboration finds an isosinglet mass of 1390(30) MeV [H from 



one recent analysis and 1380(40) MeV @ from another. Mark III finds 1430(40) MeV |6 
We take the mass of the physical mixed state with largest contribution coming from nn 
to be 1404(24) MeV, the weighted average of 1390(30) MeV and 1430(40) MeV. The mass 
of the physical mixed states with the largest contributions from ss we take as the mass of 
/o(1500), for which the Particle Data Group's averaged value is 1505(9) MeV. The mass 
of the physical mixed state with the largest contributions from the glueball we take as the 
Particle Data group's averaged mass of / (1710), 1697(4) MeV. 

Adjusting the parameters in the matrix to give the physical eigenvalues we just specified, 
m g becomes 1622(29) MeV, m S s becomes 1514(11) MeV, and E(m s ) becomes 64(13) MeV, 
with error bars including the uncertainties in the four input physical masses. The unmixed 
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trig is in good agreement with world average valence approximation glueball mass 1632(49) 
MeV, and E(m s ) is consistent with our measured value of 43(31) MeV. 
For the three physical eigenvectors we obtain 

|/ (1710) >= 

0.859(54)|^ > +0.302(52)|ss > +0.413(87)|rm >, 
|/o(1500) >= 

-O.128(52)|0 > +0.908(37)|ss > -0.399(113) \nn >, 
|/o(1390) >= 

-0.495(118)1^ > +0.290(91)|ss > +0.819(89) \nn > . 

The mixed / (1710) has a glueball content of 73.8(9.5)%, the mixed / (1500) has a glueball 
content of 1.6(1.4)% and the mixed /o(1390) has a glueball content of 24.5(10.7)%. Since, 
as well known, the partial width T(J/^> — > 7 + h) is a measure of the size of the gluon 
component in the wave function of hadron h, our results imply that T( J/\l/ — > 7 + /o(1710)) 
should be significantly larger than r(J/* -> 7 + / (1390)) and T(J/V -> 7 + / (1390)) 
should be significantly larger than T(J/ty — > 7 + / (1500)). These predictions are supported 
by a recent reanalysis of Mark III data ||16|| . In addition, in the state vector for / (1500), 
the relative negative sign between the ss and nn components will lead, by interference, to 
a suppression of the partial width for this state to decay to KK. Assuming SU(3) flavor 
symmetry for the two pseudoscalar decay coupling of the scalar quarkonium states, the 
total KK rate for / (1500) is suppressed by a factor of 0.39(16) in comparison to the KK 
rate for an unmixed ss state. This suppression is consistent, within uncertainties with the 
experimentally observed suppression. 
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FIGURES 




FIG. 1. Scalar quarkonium mass as a function of quark mass for (3 of 5.93. 
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FIG. 2. Lattice spacing dependence and continuum limit of the scalar ss mass, continuum 
limit of the scalar glueball mass, and one sigma upper and lower bounds on observed masses. 
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FIG. 3. Glueball-quarkonium mixing energy as a function of quark mass for (5 of 5.93. 
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FIG. 4. Lattice spacing dependence and continuum limit of the glueball-quarkonium mixing 
energy E(m s ) and of the ratio E{m n )/E{m s ). 
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